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ABSTRACT 


" Model  tosto  havo  been  performed  to  determine  the  blast 

♦ 

wave  propagation  in  the  tunnel  system  of  underground 
ammunition  storage  sites  v/ith  connected  storage  chambers. 
A total  of  10  different  configurationo  were  tooted  with 
variable  chamber  volumes  (300  - 15200  anglco  be- 

tween the  branch-  and  main  paooage-way  (35*  - 90*),  and 
ration  between  tho  crooa  occtiono  of  the  branch-  and 
main  passage-way  (0,125  - 0,5).  The  nodclc  were  in 
linear  scalds  1:40  to  1:100  of  typical  full  scale  in- 
stallations. By  analysing  the  data,  it  was  possible  to 
obtain  relatively  simple  scaling  relationships  which 
contain  all  of  the  moot  important  geometrical  parameters. 


The  model  data  are  compared  with  ono  large  scale  test 
and  fair  agreement  is  found. 


\ 


y 
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1 . INTRODUCTION 

This  report  is  the  fourth  in  a series  of  five  / 1 — 5 / 
describing  the  results  from  an  extensive  series  of 
model  tests  on  underground  ammunition  storage  in  case 
of  accidental  explosions.  In  two  preceeding  reports 
/ 2,  3/,  measurements  were  presented  on  the  pressures 
in  the  storage  chamber  and  the  air  blast  in  the  tunnel 
system  of  a single  chamber  storage  site,  respectively. 
The  present  report  deals  with  the  problem  of  blast 
wave  propagation  in  the  tunnel  systems  of  connected 
chamber  storage  sites. 

The  prevention  of  propagation  of  an  explosion  from  one 
chamber  to  another  is  of  major  concern  in  the  planning 
of  connected  chamber  storage  sites,  and  protective 
blast  doors,  blast  walls,  etc.  of  sufficient  strength 
are  required. 

The  results  from  the  present  model  tests  will  provide 
the  designer  blast  data  for  engineering  use  for  a 
wide  range  of  typical  storage  site  configurations. 

The  present  results  will  also  serve  as  a starting  point 
in  the  determination  of  external  safety  distances  as 
these  in  principle  can  be  determined  from  the  blast  at 
the  exit  of  the  main  passageway.  This  problem  will  be 

discussed  in  a later  report. 

# 

The  basic  scaling  laws  used  in  the  present  work  will 
be  discussed  in  Sec.  2,  followed  by  a short  presenta- 
tion of  the  experimental  details  in  Sec.  3.  The 
experimental  results  will  be  presented  and  discussed  in 
Sec.  4,  and  the  principal  results  are  finally  summa- 
rized in  Sec.  5. 
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2.  THEORETICAL  CONSIDERATIONS 

2. 1 Scaling  Relationships 

The  basis  for  the  scaling  relationships  proposed  in 
the  present  report  is  the  familiar  Hopkinson’s  scaling 
laws  /6/,  which  were  reviewed  in  Report  I /I/.  In  the 
analysis  of  the  results  for  single  chamber  storage 
sites  in  Report  III  A /3/,  various  geometrical  parameters 
were  incorporated  in  a nondimens ional  form  and  provided 
quite  general  scaling  relationships.  This  technique 
will  also  be  attempted  for  the  present  experiments. 

Fig.  2.1  shows  a cross-section  of  a typical  connected 
chamber  storage  site  including  also  the  definitions  of 
the  parameters  believed  to  be  of  primary  importance. 
Simple  similtude  analysis  coupled  with  empirical  observa- 
tions of  the  present  test  results,  have  produced  the 
following  scaling  relationships  for  peak  pressure  p, 
impulse  I,  and  positive  duration  t+: 


p = f(QX+/Vt>  A.’/Ak) 

IAk/Q  = g(QX+/Vt,  A.’/Ak) 

t+Ak/Q  = h (QX+/Vt>  A.’/Ak) 


(2.1a) 


(2.1b) 


(2.1c) 


Referring  to  Fig.  2.1, 


(2. Id) 


or  combined  volume  of  the  storage  chamber,  V., 

f t 1 


! t 


the  branch  passage-way,  A.  L , and  twice  the 
volume  of  the  main  passage-way  out  to  the  observa- 
tion point,  Ak2L.  Furthermore,  X+  expresses 


the  energy  distribution  in  the  two "directions  of 
the  main  passage-way. 
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Earlier  tests  111  have  indicated  that 
X+  = 2(1  - vm/l80°) 

X_  = V90°  (2.1e) 


X+  is  simply  a short  form  to  express  that  either 
X”  or  X_  is  to  be  used  for  the  evaluation  of  the 
blast  wave  parameters  in  the  positive  or  the  nega- 
tive direction  of  the  tunnel,  respectively  (see  def. 
in  Fig.  2.1).  It  can  be  seen  that  the  definitions 
in  Eq.  (2.1e)  satisfy  the  boundary  conditions  in 
that 

X+  = X_  for  vffi  = 90°, 

i.e.  symmetric  energy  distribution,  and  for 

v = 180°: 
m 


X = 1 


for  vm  = 0° 
m 

1 

0 

♦ 

two  particular  cases,  Eqs.  (2.1a)  - (2.1c) 
reduce  to  the  scaling  relationships  employed  in 
Report  III  A for  single  chamber  storage  sites. 

Thus,  QX+/V^.  in  the  proposed  scaling  relationships 
represents  an  effective  loading  density  including 
not  only  the  chamber  volume  and  tunnel  volumes,  but 
also  the  energy  distribution  for  variable  angles  v . 

Obviously,  the  cross-sectional  areas  of  the  branch- 
and  main  passage-way  wil  also  effect  the  blast  propa- 
gation and  this  is  represented  by  the  ratio,  A.  /A  . 

*J  c 


Likewise, 

X+  = 
X_  = 

For  these 
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The  functional  form  of  the  proposed  scaling  rela- 
tionships has  to  he  established  by  the  experimental 
results,  which  will  be  discussed  in  Sec.  4. 

2.2  Non-Scaling  Energy  Losses 

The  simple  scaling  laws  proposed  in  Sec.  2.1  do  no 
account  for  energy  dissipating  effects  such  as  ther- 
mal energy  transmission  to  the  walls,  elastic  or 
inelastic  deformation  of  the  walls,  and  viscous 
loss  due  to  wall  roughness.  Of  these  effects,  only 
the  loss  due  to  wall  roughness  is  espected  to  in- 
fluence significantly  the  blast  propagation,  but 
then  only  over  long  distances  and/or  on  the  case  of 
large  roughness.  This  is  discussed  in  Report  I. 

According  to  a model  originally  proposed  by  Porzel 
/8/,  the  peak  pressure  attenuation  may  be  expressed 
as: 


Y(p)  = constant  - 2e(e/D)(L/D)  (2.2a) 

Here,  Y(p)  is  an  impedance  function  defined  as 

Y(p)  = ^ lnp  - ^ In  (p  + 7)  - ~ (2.2b) 

with  p the  pressure  at  a distance  L/D  along  the 
tunnel.  The  ratio  e/D  is  the  average  relative 
wall  roughness  and  an  efficiency  factor  in  the 
range  l/2^e^  1. 

• For  the  relatively  smooth-walled  steel  tubes  used 
in  the  present  tests,  the • effects  of  the  wall  rough- 
ness is  expected  to  be  relatively  small,  and  will 
be  discussed  in  more  detail  in  Sec.  4. 

3.  EXPERIMENTAL  DETAILS 

Details  of  the  models,  instrumentation,  and  data  re- 
duction are  to  be  found  in  Report  I /l / , and  we 
shall  therefore  only  summarize  the  main  feature  here. 
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The  test  programme  was  designed  to  determine  the  blast 
wave  propagation  in  the  main  passageway  of  typical 
connected  chamber  storage  sites.  A total  of  nine 
combinations  of  chamber  volumes  and  tunnel  cross 
sections  were  used  as  shown  in  Fig.  3. a.  Most  load- 

*Z 

ing  densities  varied  between  3 and  70  kg/nr  and 
pressure-time  history  was  recorded  in  most  cases  at 
six  different  distances  from  the  storage  chamber, 
ranging  from  approximately  2 to  80  tunnel  diameters. 

The  TNT  charges  used  in  the  tests  were  suspended  in 
the  middle  of  the  detonation  chamber  and  initiated 
with  electrical  blasting  cap  no.  8.  This  has  an 
equivalent  TNT  weight  of  1,5  - 0,5  g.  The  smallest 
charges  used  in  the  tests  were  8 g TNT  and  the  un- 
certainty in  the  loading  density  stemming  from  the 
blasting  cap  is  therefore  less  than  6$. 

To  measure  the  pressure-time  history  in  the  models, 
standard  measurement  techniques  were  used  and  the 
blast  wave  parameters  were  evaluated  using  special 
computer  programmes.  The  total  uncertainty  in  the 
peak  pressure  was  estimated  to  be  about  10$,  whereas 
the  uncertainty  in  the  impulse  was  estimated  to  be 
about  30$. 

4.  EXPERIMENTAL  RESULTS  AND  ANALYSIS 

A total  of  approximately  1 500  blast  wave  parameter 
values  (peak  pressure,  impulse  and  positive  dura- 
tion) were  obtained  for  the  18  configurations  dis- 
cussed in  Sec.  3.  The  results  from  tests  reported 
earlier  on  connected  chamber  storage  sites  111  have 
also  been  included  in  the  present  report.  These 
earlier  data  combined  w ith  the  present  results  thus 
provided  the  18  configurations  in  Fig.  3. a.  Only 
the  most  significant  results  will  be  presented  here 
since  part  B of  this  report  contains  all  the  pressure- 
time recordings  and  tabulated  blast  wave  parameters 
except  for  pressure-time  recordings  from  the  earlier 
tests  111 . 
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The  purpose  of  the  analysis  that  follows  are  (1)  to 
determine  qualitatively  and  quantitatively  the 
dependence  of  peak  pressure,  impulse  and  positive 
duration  according  to  the  scaling  relationship  pro- 
posed in  Sec.  2,  and  (2)  to  establish  empirical  rela- 
tionships that  identify  the  more  important  factors 
taht  can  be  applied  to  blast  predictions  for  under- 
ground storage  sites  geometrically  similar  to  the 
present  range  of  models.  It  should  be  emphasized 
that  the  final  analysis  are  empirical  and  statistical 
in  approach. 

4.1  Peak  Pressure  Scaling  Relationships 

The  basic  quantity  of  interest  in  the  present  model 
tests  is  the  peak  side-on  overpressure  p in  the 
passageway.  In  Fig.  4.1a,  p is  plotted  versus  L/D 
with  the  loading  density  Q/V^  as  parameter  for  one 
typical  configuration.  Fig.  4.1b  shows  the  same  data 
plotted  versus  Q/V^  with  L/D  as  parameter.  These  plots 
are  in  socalled  dimensionless  form,  i.e.  independent 
the  linear  scaling  factor  n.  (See  Report  I /I / ) . 

Letting  subscripts  F and  N denote  full  scale  and 
model  respectively, 

Pp  = * (4.1a) 

and  Lp/Dj,  = nL^/nD^  = L^/b^ 

For  the  loading  density, 

Qp^F  = n = (4.1b) 

The  results  in  Figs.  4.1a  and  4.1b  may  therefore 
in  principle  be  used  for  any  geometrically  similar 
systems.  However,  it  is  of  interest  to  reduce  the 
number  of  curves  necessary  to  characterize  the 
blast  wave  propagation  in  this  particular  system. 

In  Sec.  2.1,  it  was  proposed  to  scale  the  data 
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according  to  QX  /V.  , where  now  V^.  is  the  total 

T U O 

volume  of  the  chamber  and  the  tunnel  out  to  the 
observation  point.  Pig.  4.1c  shows  the  results 
in  this  presentation  for  the  same  data  as  before, 
notably  with  X+  = 1 according  to  the  definition  in 
Eq.  (2.1e).  As  may  be  seen,  the  data  points  are 
evenly  distributed  around  one  straight  line.  This 
means  that  one  universal  curve  appears  to  reproduce 
the  tunnel  pressure  at  any  point  in  the  tunnel  for 
any  charge  Q for  the  present  range  of  measurements. 

Pigs.  4. Id  - 4.1i  show  the  variation  of  p versus 
QX  /V.  for  all  the  18  configurations  in  the  tests 

+ U t 

wiih  v , A.  /A,  , and  V.  as  parameters.  In  Fig.  4.1j, 

IQ  J iC  1 

the  peak  pressure  data  are  compared  for  four  widely 
different  chamber  volumes. 

Prom  Pigs  4. Id  - 4. If  it  may  be  seen  that  there  is 

a systematic  decrease  in  pressure  for  decreasing  • 

values  of  A.  /Av.  The  results  in  Figs.  4.1g  - 4.1i 
J * 

show  that  p is  on  the  average  a unique  function  of 
QX+/V.j.  for  different  values  of  the  angle  v . 

In  Fig.  4.1;),  the  peak  pressure  data  are  compared 
for  four  widely  different  chamber  volumes  V.  with 

1 . . Q 

Aj  /A^  = 1/2  and  vffl  = 90  . As  may  be  seen,  the 
data  are  relatively  insensitive  to  the  value  of  V. . 

The  qualitative  observations  made  here  will  be 
tested  more  systematically  in  the  following  sections 
which  contains  the  results  from  various  empirical 
fits. 


4.1.1  Comparisons  of  Geometrically  Similar  Models 

Of  the  18  different  models  listed  in  Pig.  3. a, 
only  two  were  exactly  geometrically  similar: 


Model  1 


V.  = 800  cm5 

i o 

A.  = 10  cm 

J 

A^  = 20  cm2 


m 


90° 


9 


Model  2 


vi 


I 


m 


15200  cm5 
70  cm2 

140  cm2 

90° 


The  linear  scaling  factor  between  Model  1 and  2 is 

n = 2,6  and  in  linear  scales  1:100  and  1:38  of 

2 

typical  full  scale  sites  with  = 20  m . 

Figs.  4.1.1a  and  4.1.1b  compares  the  peak  pressure 
data  at  two  scaled  distances  L/D  =11  and  L/D  = 20, 
respectively.  There  is  good  agreement  between  the 
two  sets  of  data  in  Fig.  4.1.1a,  but  in  Fig. 4. 1.1b 
the  pressures  for  Model  1 appear  to  be  consistently 
lower  than  those  for  Model  2.  This  difference  is  in 
fact  consistent  with  the  effects  of  wall  roughness 
attenuation  according  to  the  model  discussed  in 
Sec.  2.2.  Letting  pQ  denote  the  fictitious  "smooth- 
walled"  pressures  and  p^  and  p2  the  recorded  pressures 
for  Model  1 and  Model  2 respectively,  Eq.  (2.2a)  can 
be  expressed  as: 


Y(P0)  ~ Y(Pl)  = 2 6 e/D1 (L/D)  (4.1.1a) 

for  Model  1 (D^  = 50  mm),  and 

Y(p0)  - T(p2)  = 2€?/D2(L/D)  (4.1.1b) 

for  Model  2 (l>2  = 134,5  mm). 

The  effective  average  wall  roughness  2 £ e in 
units  of  mm  is  assumed  to  be  the  same  in  both 
models . 

Thus,  from  Eqs  (4.1.1a)  and  (4.1.1b), 

Y(P2)  - Y(p^  ) = 2 6 ¥(  1 /D.|  - 1/D2)  (L/D)  (4.1.1c) 

This  coupled  with  the  results  in  Fig.  4.1.1b 
produces 


10 


26  e ^ 0,8  - 0,5  mm 

As  € a 1 for  small  L/D- values  / 1 / , e n 0,4  - 0,3  am., 
which  is  consistent  with  the  direct  measured  wall 
roughness. 

These  limited  sets  of  data  thus  rule  out  the  presence 
of  significant  non-scaling  energy  losses  discussed 
in  Sec.  2.2.  (conduction,  radiation,  elastic  or  in- 
elastic deformation  of  the  walls). 


4.1.2  Test  of  Energy  Distribution 

It  was  postulated  in  Sec.  2.1,  Eq.  (2.1a),  that  the 
peak  pressure  will  he  a unique  function  of  the  effective 
loading  density,  QX  /V . , which  contains  the  energy 
distribution  factor  "X  defined  in  Eq.  (2.1e). 


As  an  example,  Fig.  4.1.2a  shows  the  data  for  one 
particular  case  with  vffi  = 35°  to  check  the  adequate- 
ness of  this  proposed  energy  division.  Under  the 
assumption  of  a linear  functional  dependence  in  a 
log-log  plot 

C? 

P = C1(QX+/Vt)  % (4.1.2a) 


least  squares  fits  of  the  two  sets  of  data  pro- 
duced the  parameters: 


Cx  = 7,1  - 1,1 
C9  = 0; 52  - 0,05 


for  X + 


C1  = 6,2  - 1,1 
C2  = 0,53  1 0,02 


J 


for  X - 


Clearly,  for  these  particular  sets  of  data,  the 
proposal  energy  distribution  is  adequate.  Similar 
checks  for  the  remaining  5 configurations  in  Fig.  3. a 
with  vffi  = 90°  provided  similar  results. 


4.1.3  Empirical  Scaling  Relationships 

This  section  lists  three  methods  examined  for  scaling 
the  peak  pressure  data. 
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Method  a 

Prom  the  empirical  observations  exemplified  in  the 
preceeding  sections,  it  is  postulated  that  the  peak 
pressure  may  be  expressed  as  (A.  /A,  = constant): 

Bp 

p = B1(QX+/Vt)  ^ (4.1.3a) 

where : 

B.j  = scaling  factor 
Bg  = scaling  exponent 

Least  squares  fits  of  Eq  (4.1.3a)  were  made  to 

the  data  for  each  of  the  18  configurations  in 

Pig.  3. a.  Pigs.  4. Id  - 4 . 1 i display  all  the  data 

together  with  curves  resulting  from  this  fitting 

procedure.  The  corresponding  scaling  parameters 

B.j  and  B2  are  given  in  Table  4.1.3a,  from  which  it 

can  be  seen  that  there  appears  to  be  no  systematic 

variations  in  B.  and  B9  for  different  chamber 

' 1 . 

volumes  for  the  same  value  of  A.  /Av.  An  exception 
is  the  results  for  the  = 300  cm3  chamber  volume, 
which  produced  significantly  higher  B^ -values  and 
lower  B2~values  than  those  for  the  other  configura- 
tions. It  may  also  be  noted  from  Table  4.1.3a  that 

t . 

for  decreasing  values  of  A^  /Aj,,  there  is  on  the 
average  a distinct  systematic  decrease  in  the  value 
for  B^ , whereas  B2  appears  to  show  little  variation. 

On  the  whole,  the  empirical  fits  reproduce  the  data 
remarkably  well  as  seen  from  the  error  column  in 
Table  4.1.3a.  The  average  error  for  all  measurements 
is  approximately  14$.  This  is  quite  satisfactory 
considering  the  estimated  10$  experimental  uncertainty 
in  the  peak  pressure  determination. 

Method  b 

The  second  method  used  to  scale  the  peak  pressures 
included  the  effects  from  the  variation  of  A^  /A^: 

C C 

p - c^qx/it*)  2 (Aj’/V  3 


(4.1.3b) 
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This  equation  was  fitted  to  all  18  configurations 
with  , C 2 , and  as  adjustable  parameters.  The 
results  of  this  procedure  are  shown  in  Table  4.1.3b. 

Method  c 

The  third  method  used  to  3cale  the  peak  pressures 

included  also  the  effects  from  the  variation  of 

A.  /A,  and  possible  effects  from  wall  roughness 
J K 

attenuation: 


D, 


D, 


P0  = V^V  (A.  /Ak) 
Y(p0)  - Y(p)  = (D4/D)(L/D) 


(4.1.3c) 

(4.1.3d) 


where : 


p = fictitious  peak  pressure  without 
wall  friction 

p = measured  peak  pressure 

Y(pQ)  = impedance  function  of  pQ,  Eq.(2.2b) 

Y(p)  = impedance  function  of  p,  Eq.  (2.2b) 

, I>2,  = scaling  parameters 

D,  = 2 £ e = effective  wall  roughness,  see 
4 Report  I. 


The  data  for  all  18  configurations  were  fitted 
to  Eqs.  (4.1. *30)  and  (4.1.3d).  The  results  of  this 
procedure  are  shown  in  Table  4.1.3c.  As  may  be  seen, 
the  value  for  is  rather  small,  but  negative,  which 
is  not  consistent  with  the  proposed  model.  However, 
within  the  stated  error  limits,  may  in  fact  be 
zero,  and  the  fit  therefore  shows  that  the  average 
wall  roughness  must  be  small  in  accordance  with  the 
use  of  the  relatively  smooth-walled  steel  tubes  in 
tests. 

4.1.4  Evaluation  of  the  Empirical  Fits 

On  the  basis  of  the  results  from  the  three  different 
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empirical  fits  found  in  the  preceeding  section, 

there  are  three  major  conclusions  which  may  be  reached: 

i)  Induvidual  fits  to  each  of  the  18  configurations 
(Method  a)  reproduce  the  data  within  approximately 
the  expected  experimental  scatter  (10$). 

ii)  Method  b)  excluding  the  effects  from  wall  rough- 
ness and  method  c)  including  these  effects  do  not 
produce  significantly  different  scaling  parameters. 
This  shows  that  the  tube  wall  roughness  is  small. 

iii)  Method  a)  reproduces  the  data  better  (14$)  than 
method  b)  (20$).  However,  the  convenience  of 
being  able  to  represent  all  the  data  in  one  formula 
favours  the  use  of  method  b). 

As  a general  conclusion  it  is  therefore  recommended 
to  use  the  results  from  method  b)  in  Table  4.1.3b  in 
the  prediction  of  the  smooth  wall  peak  pressure,  i.e.: 

p = (10,0  - 0,6)  (QX  /v  )0,55±0,°2(a//a  )0,  61  -C  , Oo 

+ "0  3 c 

(4.1.4) 

Here,  a standard  deviation  (96$  confidence 
interval)  has  been  used  on  the  scaling  parameters. 

There  are  various  limitations  imposed  on  the  present 
results  for  which  the  simple  relationship  in  Sq.(4.1.4) 
is  expected  to  be  invalid.  This  has  been  discussed 
in  Report  III.  In  particular,  for  significat  tunnel 
wall  roughness,  Eq.(4.1.4)  may  be  supplemented  by 
the  use  of  the  impedance  function  as  outlined  earlier 
in  this  section  and  in  Report  I / 1 / . As  has  been 
discussed  elsewhere  /9/,  the  wall  roughness  attenua- 
tion may  be  considerable  for  typical  tunnels  blasted 
out  of  rock,  and  the  use  of  Eq.  (4.1.4)  will  be  quite 
conservative  or  on  the  safe  side. 

4 . 2 Scaling  Relationships  for  Impulse 

The  results  for  the  impulse  are  much  more  difficult 
to  combine  empirically  in  a straightforward  way  than 


14 


those  for  the  front  pressure.  This  is  partly  due 
to  the  significantly  larger  uncertainties  involved 
here  and  special  features  are  not  easily  discernible. 
This  is  clear  from  the  typical  results  in  Fig.  4.2a. 

On  the  average,  the  impulse  decreases  slowly  with 
distance  from  the  storage  chamber,  but  there  is  no 
systematically  large  drop  near  the  chamber  as  was 
observed  for  front  pressure. 

To  test  the  proposed  scaling  hypothesis  in  Sec.  2.1, 
Eq.  (2.1b),  the  results  in  Fig.  4.2a  are  shown  in 
Fig.  4.2b  as  I A^/Q  versus  Q/V^.  To  a first  approxi- 
mation the  results  fall  on  discrete  curves  with  L/D 
as  a parameter.  Part,  or  all  of  the  L/D-dependence 
may  be  due  to  the  wall  roughness,  but  so  far  it  has 
not  been  possible  to  account  for  this  theoretically. 

Figs.  4.2c  and  4. 2d  show  the  scaled  impulse  data 
for  the  two  geometrically  similar  models,  Model  1 and 
Model  2,  which  were  discussed  earlier  in  Sec.  4.1.1. 

As  may  be  seen,  there  is  good  agreement  between  the 
two  sets  of  data. 

Fig.  4.2e  shows  a test  of  the  proposed  energy  dis- 
tribution represented  by  X+.  As  may  be  seen,  at 
the  same  scaled  distance  L7D,  the  scaled  impulse 
data  appear  to  be  adequately  reproduced  by  QX+/V\  . 

The  same  conclusion  may  also  be  reached  from  the 

results  shown  in  Fig.  4.2f  for  different  values 

of  v. 
m 

In  Fig.  4.2g,  scaled  impulse  data  are  shown  for  four 
widely  different  chamber  volumes  V.  with  A.  /A,  =1/2, 

i j k 

vffi  = 90  , and  L/D  = 21,6.  As  may  be  seen,  the 
data  are  relatively  insensitive  to  . 

Finally,  Fig.  4.2h  shows  an  example  of  the  reduction 
in  scaled  impulse  with  decreasing  values  of  the 

i . 

ratio  A.  /A^. 
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As  for  front  pressure,  various  empirical  scaling 
relationships  with  Eq.  (2.1b)  as  basis  were  fitted 
to  the  experimental  results.  In  particular,  a reason- 
able fit  to  all  the  data  was  obtained  using  the 
following  expression: 

E E 

I A^Q  = E1(QX+/Vt)  2(A'd/Ak)  3exp(-E4  L/D) . . . (4.2) 

Here,  E2,  E and  E^  were  introduced  as  scal- 
ing parameters.  The  result  of  this  procedure  are 
shown  in  Table  4.2.  As  expected,  the  scaling  para- 
meters are  relatively  uncertain,  but  Eq.  (4.2)  re- 
produces all  the  data  with  an  average  scatter  of 
approximately  30$.  This  is  quite  satisfactory  con- 
sidering the  relatively  large  uncertainties  in- 
volved in  the  interpretation  of  the  impulse. 

The  simple  relationship  in  Eq.  (4.2)  is  expected 
to  be  invalid  for  significant  wall  roughness  as 
for  tunnels  blasted  out  of  rock.  The  use  of 
Eq.  (4.2)  in  this  case  will  probably  produce  re- 
sult which  are  quite  conservative  or  on  the  safe 
side . 

4.3  Scaling  Relationships  for  Positive  Duration 

The  result  and  the  interpretation  of  the  positive 
duration  data  show  similar  characteristics  as  for 
the  impulse  data  discussed  in  Sec.  4.2.  The  comments 
made  there  essentially  also  apply  here.  As  for  the 
impulse  data,  a reasonable  fit  obtained  using  an 
expression: 

F F 

t^/Q  = F1(QX+/Vt)  2 (A!/Ak)  3 exp  (-F4  L/D) . . . (4.3) 

Here,  , Fg,  F,,  and  F^  act  as  scaling  parameters. 

The  results  of  the  nonlinear  least  squares  fit  of 
Eq.  (4.3)  to  all  the  data  are  shown  in  Table  4.3. 

This  fit  reproduces  all  the  results  with  an  average 
scatter  of  approximately  30$  which  again  is  quite 
satisfactory  considering  the  large  uncertainties 
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involved  in  the  interpretation  of  the  experimental 
results. 

4.4  Comparison  With  large  Scale  Tests 

Comparison  of  the  present  model  tests  with  a limited 
number  of  large  scale  tests  with  similar  tunnel  systems, 
indicates  that  the  proposed  midified  scaling  lav/s  are 
valid  for  engineering  purposes.  As  an  example,  some 
results  from  a large  scale  test  performed  by  the  Norwegian 
Defence  Research  Establishment /l 0/ are  included.  Using 
the  notation  in  Fig.  3.b,  the  storage  site  parameters 
were:  V.  = 26  ro? , A.  =1,5  m^,  Av  = 17  m^,  and  V = 90°. 

The  TNT  charges  detonated  in  the  storage  chamber  varied 
between  56  and  1500  kg,  and  measurements  were  performed 
at  distances  l/D  — 1-3. 

In  Fig.  4.4  we  show  the  measured  peak  pressures  versus 

Q/V. . As  may  be  seen,  these  results  are  in  fair  agree- 

ment  with  the  model  data  for  A.  /A,  = 1/8.  It  should 

J ^ 

be  noted  that  the  large  scale  storage  chamber  was  relative- 
ly small  compared  to  the  cross  section  of  the  main  passage- 
way and  the  L/D-values  were  also  relatively  small  and  in 
fact  outside  the  range  of  parameters  used  in  the  model  test 

5.  SUMMARY  AND  CONCLUSIONS 

In  this  paper  an  extensive  series  of  model  tests  on 
blast  wave  propagation  in  the  tunnel  systems  of  typical 
underground  ammunition  storage  sites  have  been  briefly 
reviewed.  Several  empirical  cross  checks  of  the  data 
from  18  different  model  configurations  have  shown  that 
the  peak  pressure  results  are  consistent  with  the  pro- 
posed scaling  laws  which  combine  a wide  range  of  im- 
portant geometrical  parameters.  In  particular,  the  peak 
pressure  in  the  main  passage-way  can  be  found  using 
an  effective  loading  density  QX+/Y^.,  where  again 

represents  the  total  volume  to  the  observation  point, 
and  the  parameter  X+  includes  the  effects  from  variable 

angles  between  the  branch  and  main  passage-way. 
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The  presence  of  non-scaling  parameters  (e.g.  tunnel 
wall  roughness)  makes  it  difficult  to  present  a formal 
basis  for  the  validity  of  the  proposed  scaling.  How- 
ever, direct  comparisons  of  the  present  results  with  a 
limited  number  of  large  scale  tests  with  similar  tunnel 
systems,  have  shown  that  the  proposed  modified  scaling 
lav/s  are  valid  for  engineering  purposes. 
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Table  4.1. 

3a  Least 
peak 

squares  fits 
pressure  data 

of  Eq.  (4 
for  each 

.1.3a)  to  the 
configuration . 

V. 

1 

(cm3) 

( cm2 ) 

Aj/flk 

Angle  v^ 
(deg) 

Bi 

B2  Errorb 

* 

300 

20 

l 

90 

10 , 3-1,1 

0,40*0,03 

13 

300 

20 

1/4 

90 

7, 0*1,1 

0,40*0,05 

19 

300 

20 

1/8 

90 

4, 7-1,1 

0,40*0,04 

15 

800 

20 

1 

90 

6 >5*1,1 

0,55*0,03 

14 

800 

20 

1/4  • 

90 

4 ,8-1,1 

0,52*0,05 

19 

800 

20 

1/8 

90 

2, 8-1,1 

0,58*0,03 

15 

1750 

20 

h 

90 

5,2*1, 0 

0,59*0,02 

10 

1750 

20 

1/4 

90 

3,3-1, 1 

0,59*0,02 

10 

1750 

20 

1/8 

90 

1, 7*1,1 

0,67*0,03 

14 

7250 

140 

5 

90 

6, 9*1,1 

0,63*0,04 

12 

10900 

140 

5 

90 

7, 1*1,1 

0,41*0 ,06 

23 

15200 

140 

I 

90 

6, 4*1,1 

0,50*0,05 

' 19 

7250 

140 

5 

60 

7, 0*1,0 

0,58*0,02 

13 

10900 

140 

\ 

60 

6, 3*1,1 

0,54*0,03 

14 

15200 

140 

\ 

60 

6, 2*1,1 

0°,  57*0, 03 

13 

7250 

140 

1 

5 

* 35 

6, 7*1,0 

0,54*0,02 

12 

10900 

140 

5 

35 

5, 9*1,1 

0,58*0,03 

12 

15200 

140 

\ 

35 

6, 4*1,1 

0,54*0,02 

13 

a.  In 

units 

where  p 

is  in  bar,  Q 

in  kg  TNT 

and  V.  in 

m3. 

b.  The  average  difference  between  the  fitted  curves  and 
•experimental  data. 


Table  4.1.3b  Nonlinear  least  squares  fit  of  Eq.(4.1.3b) 


to  the  combined  peak  pressure  data  for  all 
configurations.  The  average  difference  be- 
tween the  fitted  curves  and  the  experimental 
data  is  20$. 

Parameter 

Numerical  value 

C a 
C1 

10,0  - 0,3 

C2 

0,548  - 0,012 

C3 

0,61  - 0,03 

•5 

In  units  where  p is  in  bar,  Q in  kg  TNT,  and  in  nr . 


a 


Table  4.1.3c 

Nonlinear  least  squares  fits  of  Eqs.(4.1.3c) 
and  (4.1.3d)  to  the  combined  peak  pressure 
data  for  all  configurations. 

The  average  difference  between  the  fitted 
curves  and  experimental  data  is  20$. 

Parameter 

Numerical  value 

V 

9,9  ± 0,3 

« 

D2 

0,556  - 0,011 

D5  • 

0,67  i 0,03 

D4 

-0,19  - 0,20  (mm) 

a In  units 

where  p is  in  bar,  Q in  kg  TNT,  and  in  m . 

a 


Table  4.2  Nonlinear  least  squares  fit  of  Eq.  (4.2)  to 
tbe  reduced  impulse  data,  IA^/Q,  for  all 
configurations.  The  average  difference  be- 
tween the  fitted  curves  and  experimental  data 
is  29/S. 


Parameter 

Numerical  value 

Bi a 

20,8  - 0,4 

E2 

-0,399  ± 0,013 

E3 

0,0275  - 0,0011 

E4 

0,25  ± 0,05 

a 


In  units  where  I is  in  bar 
and  Yt  in  m^. 


ms 


, Q in  kg  TNT,  A^.  in 


m 


Table  4.3  Nonlinear  least,  squares  fit  of  Eq  (4.3)  to 
the  reduced  positive  duration  data,  -t+ A^/Q 
for  all  configurations.  The  average  difference 
between  the  fitted  curves  and  experimental 
data  is  33$. 


Parameter 

Numerical  value 

pi a 

00 

4* 

1 + 
o 

F2 

-0,90  - 0,01 

F3 

0,0349  - 0,0008 

F4 

-0,48  - 0,04 

a 2 

In  units  where  t is  in  ms,  Q in  kg  TNT,  Av  in  m , 

7 T JC 

and  V^.  in  nr. 


Q ■ HIT  charge 


o Chamber  volume 

t 

L o Length  of  branch  passage-way 

t 

Aj  o Cross-sectional  area  of  branch  passage-way 
A^  » Cross-sectional  area  of  main  passage-way 
vm  o Angle  between  branch-  and  main  passage-way 
L •»  Distance  to  observation  point 


D o "}4Aj J7r'  = Diameter  of  main  passage-way 

Vt  • ♦ A^*L  + 2 A^  « Total  effective  volume 

\ • 2(1  - v /1 08°)  * Energy  distribution  factor, 

positive  direction 

* **  v /90°  • Energy  distribution  factor,  negative 

n direction* 


Pig.  2.1.  Definitions  of  parameters  used  for  the 


connected  chamber  storage  sites. 


vi 

(cm5) 

** 

(cm2) 

A » 
d 

(cm2) 

vm 

(deg) 

V1  = 300 

A1  = 20 

A1'=  2,5 

v1  = 35 

V2  = 800 

A2  = HO 

A2‘-  5 

v2  = 60 

V5  = 1760 

= 7250 

v5  = 10900 

Vg  = 15200 

a3'=  10 

V 70 

v3  = 90 

Fig. 3. a.  Model  configurations  used  in  the  connected  chamber 

storage  tests.  A total  of  18  different 

combinations  were  tested  varying  the  chamber 

volume  V. . the  angle  v between  the  branch-  and 
i’  ° m 

main  passage-way  and  the  ratio  between  the  cross 
sections  of  the  branch-  and  main  passage-way, 

A.'  /A,  . These  configurations  were: 

Vi  A1  v3  i=1»2>5,  o=1,  2,  3 

Vi  A2  A4*  vm  with  1 = 4,  5,  6,  m = 1,  2,  3 


P ( bar  ) 


Pig#  4.1a.  Peak  pressure  versus  distance  in  tunnel 
diameters  for  various  loading  densities. 


(jDq)d 


Pig*  4.1b.  Peak  pressure  versus  chamber  loading  density 
with  L/D  as  parameter* 


p (bar) 


Peak  pressure  versus  effective  loading 
density,  v/ith  L/D  as  parameter.  The  straight 
line  represents  a least  squares  fit  to  all 
the  data* 


Pig.  4.1c. 


Pig.  4. Id.  Peak  pressure  versus  effective  loading  density 
for  three  different  values  of  the  ratio  be- 
tween the  cross-sections  of  the  branch  and 
main  passage-way.  Numbers  in  parentheses 
refer  to  a typical  full  scale  site  (n  = 100). 
The  straight  lines  represent  least  squares 
fits  to  the  data. 


Q/ Vt  ( k g / m ^ ) 

Fig.  4.1e.  Peak  pressure  versus  effective  loading  density 
for  three  different  values  of  the  ratio  between 
the  cross-sections  of  the  branch  and  main 
passage-way. 

Numbers  in  parentheses  refer  to  a typical  full 
scale  site  (n  = 100).  The  straight  lines  repre- 
sent least  squares  fits  to  the  data. 


p (bar) 


for  three  different  values  of  the  ratio  between 
the  cross  sections  of  the  branch  and  main 


passage-way.  Numbers  in  parantheses  refer  to  a 
typical  full  scale  site  (n  = 100).  The  straight 
lines  represent  the  least  squares  fits  to  the  data. 
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Fig.  4.  1g.  Peak  pressure  versus  effective  loading  density 
for  three  different  values  of  the  angle  between 
the  branch  and  main  passage-way.  The  numbers 
in  paran theses  refer  to  a typical  full  scale 
site  (n  = 38).  The  traight  line  represents  a 
least  squares  fit  to  all  the  data. 


p (bar) 


Pig*  4.1h,  Peak  pressure  versus  effective  loading  density 
for  three  different  values  of  the  angle  between 
the  branch  and  main  passage-way.  The  straight 
line  represents  a least  squares  fit  to  all  the 
data. 


P (bar) 


• Qx±/Vt(kg/m3) 

Pig.  4.U.  Peak  pressure  versus  effective  loading  density 
for  three  different  values  of  the  angle  between 
the  branch  and  main  passage-way.  The  straight 
line  represents  a least  squares  fit  all  the  data. 
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Pig.  4 . 1 D . Peak  pressure  versus  effective  loading  density- 
in  multiple  chamber  storage  sites  for*  four 
different  chamber  volumes.  With  a typical  full 
scale  value  A^.P  = 20  m2  for  the  main  passage-v/ay 
cross-section,  the  corresponding  full  scale 
chamber  volumes  are  listed  under  . 
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Pig.  4.1.1a.  Comparison  of  the  peak  pressure  data  versus 

loading  density  for  two  geometrically  similar 
models  as  discussed  in  the  text. 


Pig.  4.1.1b.  Comparison  of  peak  pressure  data  versus 

. loading’  density  for  two  geometrically  similar 
models  as  discussed  in  the  text. 
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Pig*  4.1  •2a.  Test  of  energy  distribution  in  the  t wo  directions 
of  the  main  passage-way  as  discussed  in  the  text. 
The  dashed  lines  represent  least  squares  fit3  to 
the  data  for  tho  two  directions  and  tho  solid 
line  a fit  to  all  the  data. 
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Pig*  4.2a*  Impulse  versus  distance  in  tunnel  diameters 
for  various  loading  densities. 


bar-ms*m2 


103- 


Vj  =1750  cm3 
A':/Ak=1/2 
= 90° 


D) 


10V 


L/D  = 4,0  • 

11.6  o 

21.6  * 

59.7  a 

85.7  ■ 


a 


10'  - 


i'll'"!  ' i i 1 i ml 


10 


10 


0 


10' 


Q/Vj  (kg/m^) 


10 


Pig.  4.2b.  Scaled  impulse  versus  effective  loading  density 
with  L/D  as  parameter. 
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Pig.  4.2c.  Comparison  of  scaled  impulse  data  versus  loading 
density  for  two  geometrically  similar  models  as 
discussed  inn  the  text. 
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Fig.  4. 2d.  Comparison  of  scaled  impulse  data  versus  loading 
density  for  two  geometrically  similar  models  as 
discussed  in  the  text. 
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Fig.  4.2e,  Scaled  impulse  versus  effective  loading  density 
to  test  the  proposed  energy  distribution. 
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Pig,  4*2f*  Scaled  impulse  versus  effective  loading 
density  to  te3t  the  proposed  energy  dis- 
tribution* 
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Pig*  4.2g.  Scaled  impuloe  versus  effective  loading 
density  for  different  chamber  volumes. 
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Pig*  4*2h.  Scaled  impulse  versus  effective  loading 
density  for  different  ratios  of  /A^. 
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Peak  pressure  versus  effective  loading  density- 
found  in  a large  scale  test  / 1 0/  in  comparison 
with  model  data. 


